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Abstract 



■ A Hamiltonian system with a modified Henon-Heiles potential is investigated. This 

\& ', describes the motion of free test particles in vacuum gravitational pp-w&ve space- 

times with both quadratic ('homogeneous') and cubic ('non-homogeneous') terms in 
the structural function. It is shown that, for energies above a certain value, the motion 
is chaotic in the sense that the boundaries separating the basins of possible escapes 
become fractal. Similarities and differences with the standard Henon-Heiles and the 
monkey saddle systems are discussed. The box-counting dimension of the basin bound- 
i . aries is also calculated. 
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1 Introduction 



Various open Hamiltonian systems in which unbounded orbits exhibit a chaotic behaviour 
have been studied (in particular in the context of chaotic scattering), see e.g. fj]]-f| and 
elsewhere. A standard method of investigation and description of these systems is to examine 
a typical set of initial conditions and observe which of the outcome possibilities is assigned 
to each initial data. The presence of fractal basin boundaries separating the possible escapes 
in the space of initial conditions is a suitable and invariant characterization of chaos 

Among the simplest Hamiltonians of this type are those with two degrees of freedom and 
polynomial potentials, in particular the famous Henon-Heiles (HH) system || given by 

H = ±(pl+P 2 y ) + V(x,y) , (1) 

where 

V = V HH (x,y) = \{x 2 + y 2 ) + |(x 3 - 3xy 2 ) . (2) 

In the literature, however, this is usually studied at energies below the critical escape value 
E c — ~, for which motions are bounded. These are regular for lowest E but become chaotic at 
higher energies still obeying E < E c . Chaos is also present above E c , where unbounded orbits 
choose between three possible distinct escapes to infinity (see, e.g. 0). The corresponding 
fractal basin boundaries have recently been investigated in JED] . 

In this paper a system with a potential similar to the standard HH model (0) is presented 
and studied. This naturally arises in the context of test particle motion in plane-fronted 
gravitational waves with parallel rays (pp-w&ves) [11|. It is well-known that in the simplest 



(so called homogeneous) pp-wscves the geodesic motion is regular but, as shown recently 



[0-0], it is chaotic in the non- homogeneous case. It is the purpose of this paper to discuss 
in detail a combined effect of both homogeneous and non-homogeneous gravitational-wave 
contributions. In section 2 the geodesic equations of motion in pp-waye spacetimes are 
presented in a suitable form in which the relation to Hamiltonian dynamics is established. 
Specific energy manifolds are then described and classified in section 3. The onset of chaos 
investigated by numerical simulation of orbits in an interesting medium energy region is 
described in section 4. Finally, in section 5 the box-counting dimension of the fractal basin 
boundaries is calculated. 
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2 Geodesies in gravitational pp- waves with both homo- 
geneous and non-homogeneous components 



The line element of the well-known class of gravitational pp-w&ve solutions JTTJ] in vacuum 
can be written 

ds 2 = 2 dCdC - 2 dudv -(/ + /) du 2 , (3) 

where f(u, C) is an arbitrary complex function of the retarded time u and the complex spatial 
coordinate (. As shown in [JT^] , |I3| , to determine the geodesic motion in these space-times 
it is crucial to solve the complex equation 

C + (U 2 /2) | = 0, (4) 

for C( r )> where a dot denotes the derivative with respect to the normalized 'time' parameter 
r, and U = u is a constant. (The remaining functions w(r) and v(r) can be found by 
subsequent integrations.) Introducing two real coordinates by x = Tie (, y = Tm (, the 
Eq. (f|) can be considered as a consequence of the Hamiltonian ([!]) with the potential 
V(x,y,u) = {U 2 /2)Kef. 

The space-time (||) with the function / linear in ( represents the trivial flat Minkowski 
universe. The quadratic function / ~ ( 2 describes homogeneous £>£>-waves (see [jllj and refer- 



ences therein), in which the particle motions are integrable. Geodesies in non- homogeneous 
£>^>-waves given by / ~ C\ n = 3,4, 5, • ■ ■, have only recently been investigated ||12|| , |13| . 



The corresponding polynomial potential V(x, y) ~ Tie ( n is called an 'n-saddle' (a 'monkey 
saddle' for n — 3). It has been found that motion in these gravitational waves with con- 
stant 'amplitude' are unbounded with n possible escapes to \(\ — > oo, and that the basin 
boundaries separating these outcomes are fractal. For non-constant amplitudes describing 
sandwich and impulsive gravitational waves a smearing of chaos is observed |14| . 

Now, it is natural and also physically interesting to investigate the combined effect of 
quadratic and higher polynomial terms in the function /. This is the main purpose of our 
work here. The simplest situation of this type occurs when / has both quadratic and cubic 
terms, 

/(C) = aQ 2 + pC 3 , (5) 

where a and (3 are real positive constants. With this, (f|) reduces to ( + (a( + /3( 2 ) U 2 = 0, 
in which the (non-trivial) constants a, (3 and U can be removed by a simple transformation 

C(r) - (jfr) ■ (6) 



In the real coordinates x and y, the equations of motion then take the form 

x + x + x 2 — y 2 = 0, i i — y — 2xy = . (7) 
These equations can be considered to follow from the Hamiltonian ([]]) with the potential 

V m (x,y) = \{x 2 - y 2 ) + |(x 3 - 3xy 2 ) , (8) 

which may be called a 'modified Henon-Heiles potential' since this has almost the form of 
the standard Henon-Heiles potential (0) with the exception of the negative sign of the y 2 
term. However, as will be shown below, this 'small' difference results in a behaviour which is 
qualitatively different from that observed in the well-known HH system. Note that, from a 
previous investigation of a generalized HH system V g (x,y) = ^(Ax 2 + By 2 ) + ^(Cx 3 — 3Dxy 2 ) 
'with adjustable coefficients', it is known that the system (§) (corresponding to the choice 
A = C = D = 1, B = — 1 in V g ) is not mtegrable . 

3 Classification of possible orbits according to a con- 
served energy 

First, we shall investigate specific properties of the polynomial potential @ and com- 
pare them with those of the standard HH potential (§) and the monkey saddle potential 
Vms( x ,y) — f^ 3 — 3xy 2 ), in which the quadratic terms are missing. There are two equi- 
librium saddle points of the potential V m at (x,y) = (0,0) and (—1,0). This is different 
from the HH potential (0) for which there is an elliptic point at the origin (in addition to 
three saddle points), and from the monkey saddle potential which has only one unstable 
equilibrium point (at the origin). The equipotentials V m (x,y) = E = const, of the potential 
(|5p are drawn in Fig. 1. Let us now briefly comment on various possible situations shown in 
the figure: 

• For E < the energy manifold H = E consists of three disjoint components (separated 
by a 'forbidden' grey region, see Fig. 1), each of which contains a single exit to infinity. 
The similar behaviour is observed for negative energies also in the monkey saddle 
potential V ms and in the HH system. 

• For E = the two components of the energy manifold join at (0,0), the third one 
remains separated. 
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• For < E < | the energy manifold consists of two disjoint components, one containing 
a single exit to infinity (x — ► — oo), the second containing two distinct exits (y — ► 
+00, y — > —00). This situation does not occur in the monkey saddle potential, in which 
all energy manifolds with E > are connected and contain three exits to infinity. The 
energy manifolds < E < ~ for the HH system consist of four disjoint components. 
One corresponds to bounded orbits and the remaining three are unbounded, each 
containing a single exit to infinity. 

• For E — I the two components of the energy manifold join at the point (—1, 0). 

• For E > | the energy manifold is connected and contains three exits to infinity in this 
case. The topology of energy manifolds is the same for both the HH and the monkey 
saddle systems. 

The orbits in the negative energy region cannot choose between different exits, so that no 
basin boundaries are present. On the other hand, the high energy region (E > |) provides 
three possible exits to infinity and numerical simulations show that the basin boundaries 
are fractal indicating the presence of chaos. These simulations give results very similar to 
those for the monkey saddle presented in ||13|| , or for unbounded motions in the HH model 
[ |T0f| . This is not surprising as the equipotentials of both the potential (H) and the standard 
HH potential @ are well approximated by the monkey saddle equipotentials V ms = E for 
sufficiently high energies. 

Thus, the most interesting energy region to be considered here is given by the condition 
< E < g. Since the motion for E < is not chaotic while it is chaotic in the high energy 
region, chaos has to appear somewhere between E = and E — ^. This is similar to the 
standard HH model, where the low energy orbits above E = are regular while they become 
chaotic for higher energies still below E — |. Such behaviour is not present in the 'pure 
cubic' monkey saddle potential V ms in which one can clearly distinguish only two energy 
regions, negative and positive. The existence of the medium energy region is caused by the 
presence of the quadratic terms in the potentials @ and (H). In the following section we 
shall concentrate on the behaviour of orbits in this interval. 

4 The onset of chaos in the medium energy region 

For an investigation of orbits in the medium energy region < E < | by a fractal method 
(which is useful for the characterization of a sensitive dependence on initial conditions in 
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systems with unbounded orbits) it is necessary to choose a suitable typical set of initial 
data. Here we may concentrate on motion in the region with two possible exits to infinity, 
y — > ±00 (the right white region in Fig. 1 d, e, f). Let us consider orbits starting at r = 
from rest from the 'right branch' of the level curve B r (the branch of V m (x,y) = E with 
x > 0), see Fig. 2 for typical trajectories of this type. This choice of initial conditions is also 
useful for localizing periodic orbits, which play an essential role in establishing the chaotic 
behaviour in both the monkey-saddle system and the standard HH system above the critical 
energy ||. We can conveniently parametrize the starting points on B r using the polar 
coordinate angle <fi with respect to the origin (such that \<p\ < |). Each orbit is then uniquely 
characterized by its energy E and a starting angle 0. 

In fact, we can concentrate on a smaller range of E in the medium energy region. It 
is easy to prove that for energies < E < the orbits starting from B r cannot cross the 
x-axis. Indeed, the second equation of (|7|) can be rewritten as y = y(l + 2x). Since x > — | 
in the considered component of the energy manifold for E < the acceleration y has the 
same sign as y. The basin boundary between the two possible exits ('up' and 'down') will 
therefore be given by <fi = 0, as the orbits starting from B r with <fr > will always escape 
directly to y — > 00, while the orbits with <fi < will always escape to y — > —00. They never 
cross the x-axis, so the basin boundary is not fractal. The unique orbit with <ft = is an 
unstable periodic orbit Ho with trajectory along the x-axis. 

Consequently, chaos has to appear in the energy region i < E < |. Numerical sim- 
ulations show that behaviour of orbits with small becomes much more complicated for 
energies E above a certain critical value E e [55, g]. The value of E will be discussed and 
specified later. Typical orbits for E = 0.16 are drawn in Fig. 2. In addition to the unstable 
periodic orbit Ho along the x-axis, there is an unstable periodic orbit n3 with > (see the 
bold curves in Fig. 2). Another one n 2 with < (which is not shown) can be obtained by 
symmetry. 

At this point it is convenient to introduce a simple notation. As mentioned above, there 
are two possible escapes to infinity for orbits in the region being discussed. Let us assign 
them a symbol j, where j = 2 for the exit with y —>■ — 00, and j = 3 for the exit with y — » 00 
(j = 1 is reserved for the exit x — > — 00, which is present for E > |, to keep the notation 
consistent with 0] and ||14]| ). Also, let us denote by r s the time r > at which the orbit 
reaches infinity (we put r s = 00 for orbits that never escape to infinity). 

We shall now investigate the detailed structure of the orbits for the typical energy 
E = 0.16. In Fig. 3 the functions j(<p) and t 8 {4>) are plotted. (These are analogues of 
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the scattering function and the time delay, respectively.) The basin boundaries separating 
possible outcomes in the initial condition state space are clearly fractal, which is confirmed 
here by zooming in up to the fourth level. This is an invariant evidence that the behaviour 
of orbits for energy E = 0.16 is chaotic. At each level of the fractal in Fig. 3, there are three 
peaks of r s (0). The left (middle, right) peak corresponds to orbits asymptotic (up to the 
given level) to the periodic orbit H2 (IIo, II3). Observe that a fractal structure occurs in a 
region in which there are only two possible outcomes. 

Further observations analogous to those described in |14| could be made. Instead, 



we shall concentrate here on the feature which is not present in the monkey saddle system, 
namely the dependence of j(0) on the energy E. In addition to a single-parameter 0- 
dependence in Fig. 3, it is necessary to investigate the basins for the 'up' (j = 3) and 'down' 
(j = 2) exits in the two-parameter region {(E,(f>) \ ^ < E < ~,\<f>\ < (fimax = f}- The 
result for \<f>\ < 0.1 is presented in Fig. 4, in which black color denotes j '• = 3 and white 
corresponds to j = 2. As expected, the only basin boundary for E below the critical energy 
E pa 0.13 is the x-axis (0 = 0). Above the critical energy this basin boundary bifurcates 
into three 'branches', the upper (with <fi > 0) corresponding to the unstable periodic orbit 
n 3 , the middle corresponding to n and the lower corresponding to n 2 . Moreover, the basin 
boundaries above the critical energy E are fractal. For example, the function j(<f)) shown 
in the previous Fig. 3 (which is clearly fractal) can be obtained as a section of Fig. 4 for 
E = 0.16. Fig. 4 is also suitable for a better localization of the critical energy E (from 
Fig. 3 it was only clear that E < 0.16). More detailed numerical calculations (not presented 
here) led us to the value E = 0.1364 ± 0.0001. 

The chaotic part of the medium energy region (E < E < |) in Fig. 4 contains two in- 
teresting subregions which deserve special attention: energies just below E = | and energies 
just above E = E . These subregions will be discussed below. 

Fig. 5 shows an enlargement of Fig. 4 for energies close to E — |. It is obvious that 
there is another bifurcation at E\ ~ 0.165 which is somewhat similar to that which occured 
at E = Eq. The fractal structure above E\ becomes even more complex, roughly speaking 
with five peaks in t s (4>) at each fractal level, instead of three peaks shown in Fig. 3. 

However, Fig. 5 indicates more: as the energy E approaches |, at least one other bifurca- 
tion occurs. This naturally leads us to the following conjecture: There is a infinite number of 
bifurcations in the fractal basin boundary below E — I at each of which the basin boundary 
along = bifurcates into three branches. 

According to the conjecture, the density of bifurcation points grows to infinity as E — > ~, 



where the third exit (to x — > — oo) is opened. The periodic orbit Ho vanishes while a periodic 
orbit Hi in the new 'neck' of the potential is created. Above E — | there are three exits to 
infinity and three corresponding unstable periodic orbits III, n 2 , n 3 , similar to the structure 
for the monkey-saddle systemF] (see [|lj, |13|). 



It would be interesting to investigate the transition region around the first bifurcation 
point E = Eq in detail. However, it proved to be very demanding on the accuracy of 
numerical methods. Fig. 6 indicates why: typical orbits in this energy region reach infinity 
at t s ^> 1, after 'bouncing' many times between the potential walls (r s m 400 in Fig. 6, which 
is still moderate: we were able to follow number of orbits with r s > 10 4 ). The orbit shown 
in Fig. 6 started at = 0.002 with E = 0.13669 and performed 110 bounces before escaping 
to j = 2. The numerical computation is extremely sensitive to the accuracy of the method. 
Such sensitivity is present close to any fractal basin boundary shown in Fig. 4, but it is usually 
restricted to a tiny range of initial conditions. However, it is present in the whole subregion 
bordered roughly by the rectangle E < E < E + 0.0005, \<f>\ < 0.003. This region near 
the first bifurcation and similar (but smaller) regions near the higher bifurcations represent 
orbits with the maximum level of unpredictability in the system under consideration. 

5 The box-counting dimension 

Finally, we present a more complete characterization of the topology of chaotic orbits in 
the medium energy region. In discussion so far, we have restricted our analysis to specific 
(yet important) orbits starting from rest from B r . A completely general orbit can start at 
an arbitrary (not forbidden) point with given energy E and direction of the initial impulse. 
For our purposes it is crucial to investigate the character of the basin boundaries, i.e. their 
regular vs. fractal nature. All unbounded orbits starting at points (x, y) with y > (or 
y < 0) which never cross the manifold section y = has the same certain exit j = 3 (or 
j = 2, respectively). Therefore, the basin boundaries can appear only among those orbits 
which cross y = 0. In order to obtain the basin boundary portrait it is thus sufficient to 
systematically study orbits which start from the (negative) x-axis, i.e. from points P = (x < 
0,y — 0), where x G [x min (E), 0], with initial impulse p x ,p y - (A complementary discussion 
of orbits starting at x > 0, y = would be analogous.) These orbits belong to the energy 
manifold H = \{j>1, + Py) + V m (x, 0) = E. The set of initial conditions considered can thus 
be represented by a disc in the {p x ,Py) plane in which each point characterize a specific 

1 There are two other periodic orbits II4 and II5 in the monkey-saddle system, which are more difficult to 
localize. We do not discuss their analogues in our system. 
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initial velocity x = p x ,y = p y of an orbit starting at (x, 0), where x < is the solution of 
the equation x 2 + |x 3 = IE — p 2 — p 2 . The circular boundary of the disc corresponds to 
geodesies starting at x = with maximum permitted impulse p 2 x + p 2 = 2E, whereas the 
center of the disc denotes p x — = p y and x min (E) e (—1, 0). In Fig. 7 we present typical 
discs for E = 0.16 (left) and E = 0.166 (right) in which each point representing specific 
initial data has again been colored in black or white according to the exit ('up' and 'down', 
respectively) reached by the corresponding evolved orbit. Analogous discs in the monkey 
saddle potential can be found in [0, fl3|| . 

It is obvious that the discs contain fractal basin boundaries, patterns and bifurcations 
which resemble those shown in Fig. 4 and Fig. 5. Interestingly, the geometrical pictures are 
very similar despite the fact that they describe different situations. In Figs. 4 and 5 the main 
independent parameter is E, whereas E is fixed in Fig. 7 with the main parameter being x. 

The fractal structure of the basin boundaries can rigorously be characterized by the box- 
counting dimension ||, |J. A large number of randomly chosen initial conditions represented 
by points in the above discs is evolved. For each point P its neighbourhood of radius e is 
also inspected. If not all the points in the neighbourhood belong to the same basin as P 
then the point is denoted as 'uncertain'. The fraction P(e) of uncertain initial conditions 
among all chosen points (we simulated 10000 points for each e) scales as P(e) ~ e 2 ~ d , where 
the parameter d gives the box-counting dimension of the fractal basin boundaries of the 
discs shown in Fig. 7. Numerical simulation indicates that d = 1.34 ± 0.07 for E = 0.16, 
and d = 1.43 ± 0.06 for E = 0.166. In fact, the dimension grows with E from d = 1 for 
E < £ ~ 0.1364 to d -> 1.50 ± 0.05 as E -»• \. 

6 Concluding remarks 

We have discussed in detail the Hamiltonian system with the potential (§) which is a simple 
but interesting modification of the famous Henon-Heiles potential (0). Interestingly, orbits in 
the system studied correspond to the geodesic motion of free test particles in the space-times 
describing plane-fronted gravitational waves with both homogeneous and non-homogeneous 
components in the structural function (g). We have demonstrated that for energies E > 
E c fa 0.1364 the motion becomes chaotic since the basin boundaries separating two possible 
exits to infinity have a fractal structure. For E > | the third exit opens and the resulting 
chaotic behaviour has a similar structure as in the monkey saddle potential V ms described 



in JTJ, |L3|. This 'pure cubic' potential for the monkey saddle system is obtained by setting 



9 



a = in (|5]). On the other hand, if the cubic term in the function / vanishes {(3 = 0), the 
gravitational wave is homogeneous and the equations of geodesic motion become linear, i.e. 
integrable. 
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Figure Captions 

Fig. 1 The equipotentials V m (x,y) = ^(x 2 — y 2 ) + |(x 3 — 3xy 2 ) = E for nine values of 
energy E representing all typical situations. The forbidden regions are grey. 

Fig. 2 An example of orbits starting from rest from the 'right branch' of the level curve 
B r with energy E = 0.16. The initial positions on B r are parametrized by an angle 0. 
The orbits Il and n 3 are unstable periodic orbits. 

Fig. 3 Plots of the function j(0) labelling the two possible exits to infinity and the function 
T s (<f>) denoting the time at which the orbit starting at given reaches infinity (E = 0.16 
here). The enlargements up to the fourth level confirm the fractal structure of basin 
boundaries. 

Fig. 4 The basins of the two exits, 'up' (black) and 'down' (white), in the subset of initial 
conditions corresponding to orbits starting from B r . These are uniquely characterized 
by energy E and starting angle 0. 

Fig. 5 Enlargement of Fig. 4 for energies close to E — |. The basin boundaries are fractal. 

Fig. 6 Typical geodesies (0 = 0.002, E = 0.13669) 'bounce' many times between the walls 
before escaping to infinity. 

Fig. 7 The discs {p x ,P y ) for E = 0.16 and E = 0.166 show the topology of orbits in the 
corresponding energy manifold. Each point in the disc represents an orbit starting at 
specific point x < 0, y = with a given velocity. Its colour denotes one of the two 
possible exits. Again, the fractal basin boundaries indicate chaos. 
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